A two dimensional model for ferromagnetic martensites 
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(Dated: February 2, 2008) 

We consider a recently introduced 2-D square-to-rectangle martensite model that explains several 
unusual features of martensites to study ferromagnetic martensites. The strain order parameter is 
coupled to the magnetic order parameter through a 4-state clock model. Studies are carried out for 
several combinations of the ordering of the Curie temperatures of the austenite and martensite phases 
and, the martensite transformation temperature. We find that the orientation of the magnetic order 
which generally points along the short axis of the rectangular variant, changes as one crosses the twin 
or the martensite-austenite interface. The model shows the possibility of a subtle interplay between 
the growth of strain and magnetic order parameters as the temperature is decreased. In some cases, 
this leads to qualitatively different magnetization curves from those predicted by earlier mean field 
models. Further, we find that strain morphology can be substantially altered by the magnetic order. 
We have also studied the dynamic hysteresis behavior. The corresponding dissipation during the 
forward and reverse cycles has features similar to the Barkhausen's noise. 

PACS numbers: 81.30.Kf,75.30.Gw,75.60.-d 



I. INTRODUCTION 



Mixed order parameter (OP) systems often exhibit 
unanticipated material properties that can be of tech- 
nological importance. A good example is that of the 
recently studied ferromagnetic martensite such as the 
Heusler alloy. Apart from the martensite transforma- 
tion temperature, such alloys can in principle have dif- 
ferent Curie temperatures for the parent and the product 
phases that could be higher or lower than the martensite 
transformation temperature. It is believed that an inter- 
play between the associated magnetic, elastic and magne- 
toelastic contributions leads to the unique properties of 
these alloys. For example, Heusler alloys show reversible 
deformations that can be induced by application of mag- 
netic field 0,0- Ullakko et al Q were the first to report 
the field induced strain in Heusler alloy and suggest that 
the magnetic field induced strain observed in martensite 
phases arises from the twin boundary motion rather than 
magnetostriction Similar effect has been reported in 
Fe^Pds also Later much higher strains ~ 6 — 7% 
have been reported in Huesler alloys H H H HE| ■ F° r 
this reason, ferromagnetic martensites are good candi- 
dates for actuator applications. 

A large volume of experimental work has been devoted 
to study various physical pro per ties characterizing the 
magnetoelastic transition [3 El El El El 13 E3 in- 
cluding their dependence on alloying composition as the 
latter affects the transition temperatures [rH Es| . In 
view of their technological importance, several models 
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have been proposed to explain the large reversible strain 
H H H El- Apart from this, there have been other 
studies to obtain the magnetic phase diagrams for dif- 
ferent sequences of ordering of the martensite transfor- 
mation temperature and the Curie temperatures of the 
austenite and martensite phases using mean field models 

ei Ei nam mm. 

During martensitic transformation, the higher sym- 
metry phase undergoes a spontaneous deformation to a 
lower symmetry phase, for instance from a cubic parent 
phase to a tetragonal martensitic phase as in the case 
of Ni2MnGa and FerPd^. Thus, the short crystallo- 
graphic axis is always the easy direction of magnetization 
|24| and magnetic anisotropy can be quite large in some 
cases. For example, in Ni-Mn-Ga martensites, the mag- 
netic anisotropy in the martensitic phase is (10 3 J/m 3 ) 
two orders greater than the parent austenite phase plj . 
Thus, the spin ordering is expected to change across the 
twin variants of the martensite. However, there has been 
few attempts [E| to study the changes in orientation of 
the magnetic order across the variants of the martensites, 
and also across the martensite - austenite interfaces. 

The purpose of the present paper is two fold. First, 
to propose a model that deals with the martensite trans- 
formation and the magnetic transition on a similar foot- 
ing with the aim of capturing the influence of magnetic 
order on strain and vice versa. Second, to study the 
orientation of the magnetic ordering in the martensite 
variants and in the parent phase for various possible or- 
derings of the Curie temperatures relative to the marten- 
site start temperature. This exercise will be carried 
out within the context of a recently introduced model 
for square-to-rectangle martensite transformation. This 
model has been shown to capture several unusual prop- 
erties associated with athermal martensites such as the 
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power law statistics of the acoustic emission (AE) sig- 
nals, the correlated nature of AE signals when the sys- 
tem is cycled in a small range of temperatures and the 
associated shape memory effect, and the precursors effect 
HM, H3, HM. Ha E3, |3ll |33, E3|- The model is coupled to 
the magnetic order parameter through a four state clock 
model. The clock model is used as continuous spin sys- 
tem in two dimension do not support long range order. 
While the strain order parameter is governed by a Bales 
and Gooding type of equation (34[, the development of 
magnetic order is carried out using Monte-Carlo method. 
We find that the magnetic order parameter which is gen- 
erally oriented along the short axis of the rectangular 
variant changes as one crosses the martensitic twin inter- 
face or martensite-austenite interface. We also find that 
the magnetic order influences the strain and vice versa 
substantially. We also study the 'noise' associated mag- 
netic hysteresis [3^, |3(| |37j that is similar to Barkhausen 
noise. 



A. Elastic free energy 

Here we adopt a recently introduced model for square- 
to-rectangle martensite in two dimensions. While the 
elastic free energy has contributions arising from all the 
three strain components, namely the bulk dilational, de- 
viatoric and shear strains, we has used only the deviatoric 
strain as the principal order parameter |3lL I32L l33| . The 
effect of the other two strain components are included 
only in an indirect way (see below ). The strain order 
parameter e(r) is defined as [3ll l32l l33| 



dx 



dy 



where u x and u v are respectively the displacement fields 
in the x and y directions. 

The elastic free energy F e is the sum of the local free- 
energy functional Fi,(e{f)) and a nonlocal long-range 
term F^ji(e(r)) that describes the transformation in- 
duced strain-strain interaction. The scaled form of the 
local free-energy Fl is 



II. STRUCTURAL AND MAGNETIC 
TRANSFORMATIONS 



Ft = / df 



m?)) + §(Ve(r)) 2 - *(r)e(r) 



(2) 



Let Tms and Tm/ be the martensite start and fin- 
ish temperatures respectively. Defining T C A and T C M to 
be the Curie temperatures of austenite and martensite 
phases, we have six different possible ways of ordering of 
the above temperatures. They are 

a. T C M < T Ms < T C A 

b. T Ms < T r M < T A 



c. T Ms 

M 



< T r A < T r M 



d. T C M < T c < T Ms 

e. T A < T Ma < T C M 

f. T A < T C M < T Ms 

The first four cases are situations where the austenite 
phase is the first to order ferromagnetically. For the cases 
e and /, the onset of the ferromagnetic temperature is 
determined by Tm s and T C M respectively. Condition a 
is more interesting because on lowering the temperature, 
the para to ferromagnetic transition takes place in the 
austenite phase followed by a martensitic transforma- 
tion which becomes magnetic only on further lowering 
the temperature. 



where D and a are in the scaled form. In the above 
equation, /;(e(r)) is the Landau polynomial for the first 
order transformation given by 



1 /-*6 

2 V ' 



(3) 



Here, t — (T ~ T C )/(T — T c ) is the scaled temperature. 
To is the first-order transition temperature at which the 
free energy for the product and parent phases are equal, 
and T c is the temperature below which there are only two 
degenerate global minima e = ±ej\/. It is well known that 
nucleation occurs at localized defect sites in the crystal 
such as dislocations and grain boundaries |2^. This is 
modeled by stress field [3ll l33l l3^ | represented by the 
third term in Eq. [21 This modifies the free-energy /; in 
such a way that the austenite phase is locally unstable 
leading to the nucleation of the product phase. 

In principle, long-range interaction between the trans- 
formed domains results from the elimination of the shear 
and bulk strains. In our model, this term has been intro- 
duced in a phenomenological way even though one can 
use the kernel obtained from St. Venant's principle [29| . 
This has been taken to have the form 



III. THE MODEL 



F LR {e} 



drdr'G{r-r')(?(r)e z (r'). (4) 



The total free energy of a ferromagnetic martensite 
has contributions arising from the elastic energy F e (e), 
magnetic energy F m (6), and magnetoelastic coupling 
F me (e, 6) where e is the order parameter for ferro-elastic 
transformation and 9 the order parameter for the para 
to ferromagnetic transition. 



We note here that unlike the long range term emerging 
from the compatibility constraint which is linear in e, 
we have used e 2 term based on symmetry of the free 
energy under e — > —e. However, we mention here that 
we have checked that using the kernel obtained by using 
compatibility condition gives results similar to the kernel 
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used here. The kernel G(f 
Fourier representation as 



is best defined in the 



F LR {e} = \ 



dkB 



where {e 2 {r)}k is the Fourier transform of e 2 (r)- The 
quantity {e 2 (r)} fc , is the complex conjugate of {e 2 (r)} k . 
The kernel B(k/k), we use is designed to pick up the 
correct habit plane directions which in the present case 
are [11] and [11], which are the favorable directions of 
growth of the product phase. In addition, the free energy 
barrier is higher along the [10] and [01] directions. This 
is captured by the simple kernel given by 



D 



-06(k 



A)fc 2 fc 2 



(6) 



where k x and k y are unit vectors in the x and y directions. 
The step function 9(k — A) has been introduced to im- 
pose a cutoff on the range of the long-range interactions. 
The constant (3 is the strength of the interaction. The 
real space picture of B(k/k) is similar to the long-range 
interaction of Kartha et. al [29| . 

We have used Rayleigh dissipation functional (Ijj to 
account for the dissipation accompanying rapid interface 
motion described by 



(7) 



This dissipation term has been successfully used to ac- 
count for several features of acoustic energy dissipated 
during the martensitic transformation |3lL l.'S'2t l33| . 

The Lagrangian then is given by L = T — F, where 
F is the total free-energy and T is the kinetic energy 
associated with the system is given by 



T = dfp 



du x (r,t) 

at 



du y (r,t) 

at 



(8) 



Here p is the mass density. We use this for obtaining the 
equations of motion for e. 



B. Magnetic Free Energy 
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where the first term is the exchange energy, the second 
term is the magnetic anisotropic energy and the third 
the Zeeman term. The spin vector 9 is allowed to as- 
sume four values i.e 0, 90, 180 and 270° in the austenite 
phase, and four similar values with respect to the diag- 
onal in the martensite phase (see below). The sum for 
the first term is taken over the nearest neighbors. The 
coefficient j p is different for the austenite and martensite 
phase (p = a, m). k\ is the anisotropy coefficient and 
h the applied magnetic field. (<f> is the angle between 
the magnetic field and the x— axis.) In the martensite 
phase, we have used two different values for j m , one each 
corresponding to the short and long sides of the marten- 
site variant ( in addition to the anisotropy term). The 
magnetic anisotropic energy sets the orientation of the 
magnetic order parameter along the easy axis. Usually, 
one uses sin 2 9, where 9 is the angle measured with re- 
spect to the x— axis. From our earlier studies [3lIl32ll33T ]. 
we know that the martensitic variants are along the di- 
agonal. Thus, in order to have magnetic order parame- 
ter oriented along these easy directions in the martensite 
phase, we have subtracted 45°. This gives upto an addi- 
tive term sin2 9. Lastly, there is magnetostriction energy 
given by 

F roe (e(f), 9(f)) = - J dr^=e(r) sin(20(r)), (10) 

where b\ is the phenomenological magnetoelastic con- 
stant. The effect of the strain produced during the 
martensitic transformation on the magnetic transition is 
captured by this term. Here again, considering the fact 
that e takes on positive and negative values for the two 
twin variants and that these twins are oriented along the 
diagonal, the above term is minimum along the diago- 
nal. The total magnetic contribution to the free energy 
is given by the sum of F m (9(f)) (= J f m (9(f))df ) and 

± me ' 

The total free energy is minimized with respect to both 
the order parameters e and 9. Here, we note that while 
strain OP is continuous, the magnetic OP is discrete. 
Moreover, in our model, the martensitic transformation 
is considered athermal and is therefore described by de- 
terministic equations obtained using 

i = x,y. (11) 



As stated earlier, one of our main interest is to study 
the direction of magnetic order in the martensite vari- 
ants and in the parent phase. However, it is well known 
that continuous spin ( U(l) asymmetry) systems do not 
support long range order in 2-D. For this reason, we use 
a 4-state clock model. The free-energy has contributions 
from exchange interaction, Zeeman energy, and magneto- 
crystalline anisotropy energy and is given by, 



In terms of rescaled space and time variables, we get 



+ Hr,t) 
bi 



dkB(k/k){e 2 (k, t)} k e 



ik. 



v/2 



sin(26»(f)) 



(12) 



f m W)) 



Jp 
2 



^cos((9(r) - 9(r')) 



ki 



cos(20(r) - tt/4) - hcos(9{r) - 0),(9) 



Here, f3,j and b\ are scaled parameters. However, as 
the magnetic order parameter is allowed only four val- 
ues ( Ising like), it is most conveniently dealt with using 
Metropolis Monte Carlo algorithm. 
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IV. NUMERICAL SIMULATIONS 

Numerical simulations to study the evolution of strain 
morphology and the magnetic ordering are described be- 
low. We have solved Eq. ^| using Euler's scheme with 
periodic boundary conditions after discretizing it on a 
N x N grid. The grid size of the mesh used is Ax = 1 
and the smallest time step is At = 0.002. Simulations 
were carried out for N = 128 and 256. However, the re- 
sults reported here are for N — 128. A pseudo-spectral 
technique is employed to compute the long-range term 
[33l | . The cutoff A in the long-range expression is chosen 
to be 0.1. The inhomogeneous stress field a(r) is appro- 
priately chosen to describe the multi-defect configuration 
given by 




where Co (*j ) is the magnitude of the stress field at sites 
fj which are randomly chosen defect sites, jmax is the to- 
tal number of defect sites, and Q the width of the field. 
A random distribution of 0.1% defects is used which for 
N = 128 gives jmax = 16. The values of uo(fj) are taken 
to be uniformly distributed in the interval [—0.3,0.3]. 
The parameters chosen for the martensite transforma- 
tions are f3 = 30 and 7 = 1. We start our simulations 
with the system being in a homogeneous state with e(r, 0) 
distributed in the interval [-0.005,0.005]. 

The initial values of the spin Oij at the location ij are 
taken to be randomly oriented among 0, 90, 180 or 270° 
directions in the austenite phase and 45, 135, 225, or 315° 
in the martensite phase. Thus, the allowed spin flips are 
through ±90°. However, when crossing either the twin or 
the martensite-austenite boundaries, we allow spin flips 
of ±45° to ensure that the spin orientations across the 
boundaries change appropriately. The anisotropy coeffi- 
cient is taken to be k\ = 0.6 and the magnetostriction 
coefficient b\ = 0.6 in our simulations (unless otherwise 
stated). The exchange interactions j p which fixes the 
Curie temperatures of the magnetic phases is chosen so 
as to obtain the required ordering of the Curie tempera- 
tures. In the Monte-Carlo simulations, when computing 
the equilibrium magnetization at each temperature, the 
first 2.0 x 10 5 MC steps were discarded for thermalising, 
and averages are carried out over the next 1.0 x 10 5 MC 
steps. ( We find that this duration is adequate for ob- 
taining steady state values of the magnetization.) In the 
martensite phase, we have assumed that j rn is different 
along the short and long axis of the variant denoted by 
jmsh and jmio respectively. 

Before presenting our results, we note here that two 
types of simulations are possible due to separation of time 
scales associated with the two order parameters. As men- 
tioned earlier, our model describes athermal martensite, 
and thus, the time evolution of the strain order parameter 
is considered fast compared to the magnetic order param- 
eter (even though there is no direct mapping of Monte- 



Carlo time with the real time). Thus, the first possibility 
is to evolve the strain order parameter till a stationary 
morphology of the martensites is obtained and then a MC 
evolution of the magnetic order parameter is effected (or 
vice versa depending on the ordering of the transition 
and transformation temperatures). In the second type 
of simulation, the strain OP is allowed to develop for a 
short time (long enough time to reach stationary values) 
and MC evolution is carried out for about 3,00,000 steps. 
This sequence is repeated till stationary values of both 
the OP's are achieved. For most purposes, we find that 
the first type of simulations is adequate. However, we 
shall also report the results on the second type of simula- 
tion which demonstrates the substantial mutual influence 
of the two order parameters. 



V. RESULTS 

Before proceeding further, we note that it is necessary 
to make a correspondence between the temperature scale 
used in the MC simulations with the temperature scale 
used in Eq. ^] In our MC simulations, the scale of 
temperature is in the range of unity as we choose j a in the 
neighborhood of unity (in units of fcs) for the austenite 
phase. However, the scaled temperature r, for the strain 
variable ranges mostly over negative values. To map, 
we define a relation between the scaled temperature in 
MC simulations t mCi an d t such that r = maps to 
t mc = 1, ie., t mc = T ~y M ™r T T M M / ] , where t mc = T/T c 
with T c taken to be unity for a hypothetical austenite 
phase with j a = 0.885. ( tm s and tm{ are the scaled 
martensite start and finish temperature, and the value of 
TMs — TMf is determined numerically.) We have also used 
a small magnetic field along the positive x— direction to 
ensure that magnetization in the austenite phase is along 
the same direction. We now consider a few cases of the 
ordering of temperatures listed earlier. 



A. Case: T C A < T Ms < T C M 

To obtain the above ordering of temperatures, we 
choose the exchange interactions in the austenite phase to 
be j a = 0.826, and for the martensite phase jmsh = 0.95, 
jmio = 0.35 along the easy and hard directions respec- 
tively. With this, we get T C M = 1.0312 ( higher than 
Tms, ie., t ~ 0.0), which can only be determined nu- 
merically as it is influenced by b\ and k\. With the 
above value of j n , we find T C A = 0.937 corresponding 
to r = —1.98 pdjL In this case, the Curie temperature of 
the martensite phase is not detectable as T C M is higher 
than Tms and the high temperature austenite phase is 
paramagnetic. As the temperature is lowered below Tm s 
( t ~ 0.0), some regions of the austenite phase undergoes 
structural transformation to the martensite phase which 
then orders ferromagnetically. This is reflected as an 



abrupt increase in the magnetization around r = 0.0 ev- 
ident from Fig. ^ Further, magnetization grows in pro- 
portion to the area fraction of the martensite phase until 
is reached. At this point (t ps —2.0), one again finds 
a sudden increase in the magnetization(Fig. ^ •). This 
can be attributed to ordering of the untransformed frac- 
tion of the austenite phase. In addition, we note that the 
contribution to the total magnetization from the austen- 
ite phase is higher as magnetization (per unit area) of 
this phase is oriented along the x— direction compared 
to the martensite phase which is oriented along ±45°. 
However, as we decrease the temperature, the decreas- 
ing volume fraction of the austenite compensates the ori- 
ented contribution from the austenite phase to the total 
magnetization. Consequently, a further decrease in tem- 
perature leads to a marginal decrease in the magnetiza- 
tion around r = —10.0 eventually saturating to a value 
around 0.7. The sudden jump in the magnetization at 
t = —2.0 would not be seen if the value of is such that 
the volume fraction of the untransformed austenite phase 
is small. In such a case, the magnetization curve increases 
smoothly without such sharp changes. This is shown on 
the same figure ( o in Fig. ^| when Tjr = 0.765 which 
corresponds to r = —7.52. For this case, the saturation 
value is much lower than for • as the parameter values 
are slightly different (j a = 0.55, j msh = 0.90, j m i = 0.35 
and fci = 0.5). On visual inspection, we find a large num- 
ber of smaller sized twins which are not magnetized to 
the same extent as the twins for the case •. 





FIG. 2: (Online color) Morphology of the system at r = —2.0 
for the case T C M > r > T C A . Here pink and cyan corresponds 
to the two variants of the martensite and blue to the austenite 
phase. The magnetic ordering is present only in the marten- 
site variants. For the sake of visual clarity a small region of 
85 x 55 grid points is displayed. 




FIG. 1: Magnetization vs temperature plot for the condition 

Tr A _^ 'I ' rp M 



FIG. 3: (Online color) Morphology at r = —9.0 for the case 
T C M > t > T C A . The system is fully transformed to marten- 
site, and the spins are well ordered. A small region of 85 x 55 
grid points is shown. 



Figure [21 shows a morphological plot of the magnetic 
ordering in the martensite domains at r = —2.0 for the 
case shown in Fig. ^ (•). Note that the austenite phase 
is disordered. It is also clear from the Fig. [2] that the 
two different variants of the martensite phase have differ- 
ent orientations corresponding to the two easy directions 
of magnetization. Figure shows the morphology for a 
lower temperature (r = —9.0), a value at which the sys- 
tem has almost fully transformed to the martensite phase 
and the magnetization in both the twin are well ordered 
even though the twin domains are rather small. 



B. Case: T c < T Ma < T c 

As the temperature is lowered below T^, one first finds 
magnetic ordering of the austenite phase followed by the 
formation of the martensitic phase below Tm s and then 
magnetic ordering of the martensite phase. Choosing 
j a = 1.22, one can estimate to be 1.38 which cor- 
responds to = 12.28. However, in our simulations, 
the transition temperature is identified by the peak in 
the susceptibility which in this case is seen at = 12.8. 
(The shift in could arise from the finite sample size 
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0.75 



0.5 



0.25 




FIG. 4: Magnetization vs temperature plot for T c < Tms < 

rpA 
r. • 



and small but finite magnetic field value of h = 0.025 
used, ft may be mentioned here that most measurements 
of magnetization in Heusler alloys are carried on with fi- 
nite magnetic field. ) The martensite transition temper- 
ature obtained numerically is T^ 1 = 0.719 for j ms h = 0.7 
and jmio — 0.3 which corresponds to = —9.0 |41| . 



The magnetization curve for this case is shown in Fig. 
0] In our simulations, the system shows a small amount 
magnetization at r = 12.8 which corresponds to the peak 
in the susceptibility. However, one observes a rapid in- 
crease in the magnetization M as we cool further. From 
Fig. it is clear that the austenite phase reaches near 
saturation value at r = ( ~ Tms ) ■ As the system is fur- 
ther cooled below Tm s , the martensite phase is formed 
is not magnetically ordered at this temperature, which 
leads to a sudden drop in the magnetization. Indeed, 
the extent of the drop is equal to the area fraction of 
the austenite phase that transforms at r = 0. On fur- 
ther cooling below T C M , as the martensite domains start 
magnetically ordered, the magnetization again increases 
reaching a value of 0.55 as we approach r = —20. A 
snap shot of the spin ordering for r = —3.0 is shown in 
Fig. [3] It is clear that spins in the austenite phase are 
oriented along the x— axis while those in the marten- 
site phase are disordered. We note that the maximum 
value of the magnetization in the transformed marten- 
site phase (corresponding to r = —20) does not reach 
the saturation value l/v2- This can be attributed to the 
substantial area fraction in between the martensite twins 
that is not ordered. At these temperatures even as the 
martensite domains quite fragmented, spins in the two 
martensite variants are ordered along their easy axis of 
magnetization, ie., ±45°. The spin ordering here is sim- 
ilar to T C A/ < Tms < Tc & t similar temperature range 
Fig. El 




FIG. 5: (Online color) Ferromagnetic ordering of the austenite 
phase for the case T C M < T Ms < T C A at r = -3. Note that 
the spins in the austenite phase (blue color) are ordered and 
the martensite phase (pink and cyan corresponding to the two 
variants) is disordered. A small region of 85 x 55 grid points 
is shown. 



C. Case: T C A < T C M < T Ms 

To obtain the required ordering of the transition and 
transformation temperatures, we choose the exchange in- 
teraction to be j a — 0.75 and j ms h — 0.9, j m i = 0.65 for 
the austenite, and the easy and hard directions of the 
martensite phase. With this choice, we get, T C A = 0.83 
(r = -5.5 ) and T C M = 0.98 ( r = -0.64) for the austen- 
ite and the martensite phases respectively. Above T C M , 
there is no ordering. As the temperature is lowered below 
T C M , the magnetization increases till T A ( r = —5.5) at 
which one sees a small jump in the magnetization arising 
from the fraction of the untransformed austenite phase as 
shown in Fig. 6. Thereafter, the magnetization increases 
to a maximum of 0.6 and falls marginally. This decrease 




FIG. 6: Magnetization vs temperature plot for T C A < T C M < 
Tms ■ 
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can be attributed to decreasing fraction of the austenite 
phase. 



D. Influence of magnetization on strain 

It is evident from Eq. 12, that strain is affected by 
magnetization and vice versa when the system is in the 
martensite phase. This is a direct result of the magne- 
tostriction coupling. We illustrate this interplay between 
the two order parameters for the case < Tms < 
T C M . As the high temperature phase is the paramag- 
netic austenite phase, we first obtain the initial station- 
ary state for the strain order parameter by quenching 
the system below Tm s (ie., t just less than 0.0 ). In 
this simulation, keeping the initial state in the marten- 
site phase, we allow the magnetic order to develop for 
a short time ( MC steps are typically 50,000 steps) fol- 
lowed by strain development. This is repeated for several 
cycles. Both the strain and magnetic order parameters ( 
M averaged over 30,000 MC steps) are monitored at the 
end of each cycle. We find that while maximum changes 
occur in both the OP's during the first few cycles, the 
magnetic order takes longer to reach a stationary value 
going through oscillatory changes. The stationary val- 
ues of the OP's are obtained after about 12 cycles. The 
percentage change at the end of the simulation in strain 
OP, with respect to the initial value is ~ 6.0%. In con- 
trast, the magnetic order changes by almost ~ 10% with 
respect to the intial value, here taken as the value at the 
end of the first 50000 MC steps. 

As the above simulations show that there is a strong 
mutual influence of the order parameters on each other, 
we have allowed the system to stabilize in each cycle., 
ie., strain is evolved till it reaches a stationary value fol- 
lowed by the magnetic oder parameter reaching a sta- 
tionary value consistent with the strain value. ( Near 
stationary <pA values are reached in 25 time units.) Each 
magnetization run was for 3,00,000 MC steps and M is 
calculated over the last 1,00,000 MC steps. This process 
is repeated for several cycles till stationary values of the 
OP's are obtained. Again, we found that while maximum 
changes in the OP's occur in the first cycle, the mag- 
netization which goes through oscillatory changes takes 
much longer time to stabilize. The stationary values of 
the two order parameters obtained in this simulation are 
essentially the same as in the earlier simulation where the 
duration for the time development was much smaller. A 
plot of the strain and magnetization as a function of time 
and MC steps is shown in Fig. |SJ Snap shots of the strain 
morphology at the end of first and last cycle are shown 
in Fig. [7| It is evident that the martensite boundaries 
have changed significantly from the initial configuration. 
However, we note that in real systems, during time evo- 
lution, there is a continuous feed back between the two 
order parameters. In this sense, simulations where the 
OP's are allowed to evolve for short durations are closer 
to continuous time development of both the OP's. 




FIG. 7: (Online color) Strain morphology of the system for 
T C A < T M s < T C M at r = -0.5 before the start (left) of the 
cycles and at the end configuration (right). The change in the 
area fraction of the martensite phase due to the magnetization 
is clearly seen. The entire 128 x 128 grid points is shown. 
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FIG. 8: Plot of magnetization and area fraction versus MC 
steps and time for T C A < T Ms < T C M at r = -0.5. 



As the maximum change in strain occurs in the first 
cycle, we have examined the influence of temperature on 
this differential change. This is shown in Fig. [5] for two 
conditions < T Ms < T C M ( o) and < T C M < T Ms ( 
♦). As can be seen, in both cases, the differential change 
is maximum at the onset of the martensite phase. 



E. Hysteresis 

We have also studied hysteresis loop keeping the tem- 
perature at an appropriate value. Here the field is 
changed at a specified rate from zero value to a posi- 
tive maximum ( corresponding to the saturation value 
of the magnetization) and then reversing the field to the 
negative maximum, and back again to the positive maxi- 
mum completing a cycle. First, the system is thermalised 
to the desired temperature, here, r = —3.0, for the case 
T C A < T Ms < T C M . The loop is swept by varying the field 
at intervals of h — 0.05, except near saturation where it 
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FIG. 9: Difference in area fraction vs scaled temperature, 
before and after magnetization for the cases T c 4 < Tms < T^ 1 
(o) and T C A < T C M < T Ms (♦). 



FIG. 11: Dissipation plotted as a function of external applied 
field for the case T C A < T Ms < T C M at r = -3.0. (a) Reverse 
cycle (from h — 1.0 to h = —1.0 through zero filed), (b) 
Forward cycle. 



is varied at intervals of 0.2 ( beyond 0.6). Magnetization 
is evolved for 500 MC steps at each value of the field. At 
this temperature, 44% of the system is transformed to 
the martensite phase which is ferromagnetic. Figure ITUI 
shows the hysteresis cycle for this condition. 

1 i ■ ■ ■ 1 
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_1 i ' ' ' 
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FIG. 10: Hysteresis plot for the case T C A < T Ms < T C M at 
r = -3.0. 

As is clear from the previous section, magnetization 
influences strain. However, as shown in our earlier work 
[3ll l32l l33l | , strain changes are always accompanied by 
dissipation. For this reason, we have also monitored the 
dissipation accompanying strain during one cycle of the 
external magnetic field. We find that the dissipation is 
highest when the field is close to zero. Here, it is pertinent 
to recall that in our earlier work we have shown that 
the dissipation accompanying the strain changes have all 
the characteristic features of acoustic emission J^l], [32, 
l33| . If this interpretation is followed here, the dissipation 
is similar to the Barkhausen noise observed during the 
hysteresis. 



VI. SUMMARY AND CONCLUSIONS 

In summary, we have extended the recently introduced 
two dimensional square-to-rectangle martensite model to 
include the magnetic order parameter through a 4-state 
clock model to study ferromagnetic martensites. The 
model includes all the features of the structural trans- 
formation (ie., long-range strain interaction, inertia, and 
defect centers that act as nucleation sites) and magnetic 
transition (ie., the exchange interaction, anisotropy, and 
magnetostriction) . The model is used to explain the mag- 
netization curves of the austenite-martensite system un- 
der various ordering of the Curie temperatures and the 
martensite transformation temperature. By choosing ap- 
propriate exchange interactions in the two phases, we 
are able to obtain the characteristic magnetization plots. 
The model shows that the orientations of the magnetic 
order of the ±e twins are along ±45° directions, whereas 
it is orientated along the x— axis in the austenite phase. 
Further, the model also captures the abrupt change in 
the orientation of the magnetization across the two vari- 
ants of the martensite as also the austenite-martensite 
phase boundary. The model shows that there is a subtle 
interplay between the growth of the strain and magnetic 
order parameters that leads to jumps in magnetization 
at in cases when the of the austenite phase is 
lower than T C M ( or Tm s ), a feature not predicted by the 
mean field models. The model also shows that the strain 
order parameter is influenced by the magnetic order and 
vice versa in a significant way. As this mutual influence 
is triggered by the magneto-striction coupling, the model 
may have relevance to Terfenol-D. In particular, the mag- 
netic order appears to change the initial strain morpholgy 
of martensite domains in a substantial way. Studies on 
hysteresis coupled with the earlier identification of dissi- 
pation with acoustic emission shows that the dissipation 
is maximum while the applied field passes through zero 
value. This feature is similar to the Barkhausen noise. 
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The Barkhausen noise has been explained as arising from 
pinning and depinnin g of the magnetic domains. In our 
case, as shown earlier [3lL l33j . it results for the jerky 
motion of the interface. The driving force for this change 



is, however, the interaction between the strain and the 
magnetic order parameter. It is also clear that the ap- 
proach taken here can be extended to systems where fer- 
roelectric order parameter has strong coupling to strain. 



[1] A. DeSimone, and R.D. James, J. of the Mechanics and 

Physics of Solids 50, 283 (2002). 
[2] T. Yamamoto, M. Taya, Y. Sutou, Y. Liang, T. Wada, 

and L. Sorensen, Acta Mater. 52, 5083 (2004). 
[3] K. Ullakko, J.K Huang, C. Kantner, R.C.OHandley, and 

V.V. Kokorin, Appl. Phys. Lett. 69, 1966 (1996). 
[4] A. A. Likhachev and K. Ullakko, Phys. Lett. A 275, 142 

(2000). 

[5] R.D. James and M. Wuttig, Phil. Mag. A 77, 1273 

(1998) . 

[6] R.C.OHandley, J. Appl. Phys. 83,3263 (1998) 

[7] R.C. O'Handley, S.J. Murray, M. Marioni, H. Nembach, 

and S.M. Allen, J. Appl. Phys. 87, 4712 (2000). 
[8] S.J. Murray, M. Marioni, S.M. Allen, R.C. OHandley, 

and T.A. Lograsso, Appl. Phys. Lett. 77 , 886 (2000). 
[9] H. Chopra, C. Ji and V.V. Kokorin, Phys. Rev. B 61, 

R14913 (2000). 

[10] A. Steuwer, T. Mori, H. Kato, T. Wada, J. of Appl. Phys. 

94, 2761 (2003). 
[11] T. Wada, Y. Liang, H. Kato, T. Tagawa, M. Taya, and 

T. Mori, Mater. Sci. and Engg. A361, 75 (2003). 
[12] Y. Murakami, D. Shindo, M. Suzuki, M. Ohtsuka and K. 

Itagaki , Acta Mater. 51, 485 (2003). 
[13] E.C. Oliver, T. Mori, M.R. Daymond, P.J. Withers, Acta 

Mater. 51, 6453 (2003). 
[14] O. Heczko, K. Jurek, and K. Ullakko, J. Mag. Mag. Mat. 

226-230, 996 (2001). 
[15] J. Cui, T.W. Shield, and R.D. James, Acta Mater. 52, 

35 (2004). 

[16] V.A. Chernenko, C. Segui, E. Cesari, J. Pons, and V.V. 
Kokorin, Phys. Rev. B 57, 2659 (1998). 

[17] A.N. Vasil'ev, A.D. Bozhko, V.V. Khovailo, I.E. Dik- 
shtein, V.G. Shavrov, V.D. Buchelnikov, M. Matsumoto, 
S. Suzuki, T. Takagi, and J. Tani, Phys. Rev B 59, 1113 

(1999) . 

[18] V.A. Chernenko, V.A. L'vov, S.P. Zagorodnyuk, and T. 

Takagi, Phys. Rev. B 67, 064407 (2003). 
[19] D. I. Paul, J. Marquiss and D. Quattrochi, J. Appl. Phys. 

93, 4561 (2003). 
[20] V.D. Buchel'nikov, AT. Zayak, A.N. Vasil'ev, D.L. Dali- 

dovich, V.G. Shavrov, T. Takagi, and V.V. Khovailo, 

JETP 92, 1010 (2001). 
[21] V.D. Buchel'nikov, A.N. Vasil'ev, AT. Zayak, and P. 

Entel, JETP 92, 1019 (2001). 
[22] V.D. Buchel'nikov, V.S. Romanov, and AT. Zayak, J. of 

Magn. Magn. Mater. 191, 203 (1999). 



[23] A.N. Vasil'ev, V.D. Buchel'nikov, T. Takagi, V.V. Kho- 
vailo, and E.I. Estrin, cond-mat/0311433| (2003). 

[24] L. Straka and O. Heczko, J. of Appl. Phys. 93, 8636 
(2003). 

[25] R. Tickle and R.D. James, J. Mag. Mag. Mater. 195, 627 
(1997). 

[26] E. Vives, J. Ortin, L. Manosa, I. Rafols, R. Perez- Ma- 
grane and A. Planes, Phys. Rev. Lett. 72, 1694 (1994). 

[27] I.S. Muto, S. Takeda, R. Oshima, and F.E. Fujita, J. 
Phys.: Condens. Matter 1, 9971 (1989). 

[28] I.M. Robertson and CM. Wayman, Phil. Mag. A 48, 421 
(1983); ibid 48, 443 (1983); ibid 48, 629 (1983). 

[29] S. Kartha, J.A. Krumhansl, J. P. Sethna, and L.K. Wick- 
ham, Phys. Rev. B 52, 803 (1995). 

[30] R. Oshima, M. Sugiyama, and F.E. Fujita, Metall. Trans. 
A 19 A, 803 (1998). 

[31] R. Ahluwalia and G. Ananthakrishna, Phys. Rev. Lett. 
86, 4076 (2001) and references therein. 

[32] S. Sreekala and G. Ananthakrishna, Phys. Rev. Lett. 90, 
135501 (2003). 

[33] S. Sreekala, R. Ahluwalia, and G. Ananthakrishna, Phy. 

Rev. B 70, 224105 (2004). 
[34] G.S. Bales and R.J. Gooding, Phys. Rev. Lett. 67, 3412 

(1991). 

[35] E. Obrado, E. Vives, and A. Planes, Phys. Rev. B 59, 
13901 (1999). 

[36] M. Rao, H.R. Krishnamurthy, and R. Pandit, Phys. Rev. 

B 42, 856 (1990). 
[37] J. Garcia-Otero, M. Porto, J. Rivas, and A. Bunde, J. 

Mag. Mag. Mat. 203, 268 (1999). 
[38] P.L. Ferraglio and K. Mukherjee, Acta Metall. 22, 

835(1974). 

[39] Y. Wang and A.G. Khachaturyan, Acta Mater. 49, 1165 
(2001); Y.M. Jin, A. Artemev, and A.G. Khachaturyan, 
Acta Mater. 49, 2301 (2001); A. Artemev, Y. Wang, and 
A.G. Khachaturyan, Acta Mater. 48, 2503 (2000). 

[40] L.D. Landau and E.M. Lifschifz, Theory of Elasticity, 3rd 
ed. (Pergamon, Oxford, 1986). 

[41] We stress that these unrealistic values of j m sh, jmio and 
j a are necessary to have a good separation between 
and T^ 1 . The anisotropy and the magnetoelastic constant 
were taken to be ki — 0.5, and &i = 0.6 are also unre- 
alistic. These were however, the minimum values of the 
parameters for which the martensitic domains of reason- 
ably small twins sizes would get ordered along the easy 
direction in a reasonable time. 



